Quantum computation via translation-invariant operations on a chain of qubits 
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A scheme of universal quantum computation on a chain of qubits is described that does not 
require local control. All the required operations, an Ising-type interaction and spatially uniform 
simultaneous one-qubit gates, are translation-invariant. 
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I. INTRODUCTION 

Symmetry reduces complexity. In physical systems re- 
alizing quantum computers, the highest degree of sym- 
metry is therefore not the most desirable. A quantum 
computer needs to be sufficiently simple and robust to 
be controllable in an experiment yet complex enough to 
be universal. One may therefore ask the question "How 
much symmetry does a quantum computer allow for?" . 
In fact, a number of physical systems considered for the 
realization of a quantum computer such as optical lattices 
or arrays of micro-lenses [2] are translation invariant, 
and the above question acquires a practical aspect. 

Quite surprisingly, it turns out that universal quan- 
tum computation can tolerate a fair amount of symme- 
try. Recently, a scheme of quantum computation using 
the rotation-invariant measurement of the total 'spin' of 
two qubits as the only gate has been devised 0] • Further- 
more, a translation-invariant computation scheme has 
been described p|. 

The computational power of translation-invariant or 
nearly translation invariant quantum systems was re- 
vealed in Lloyd's |j| and Watrous' @ work on quantum 
cellular automata (QCA). In [|| it was shown that a one- 
dimensional QCA can simulate any quantum Turing ma- 
chine. Translation invariance is broken only by the ini- 
tial state which encodes the program. The schemes [j| 
and allow to simulate quantum circuits using a chain 
of qubits with a generic translation-invariant interaction. 
They require different species of qubits in a periodic ar- 
rangement and local addressability at one end-point of 
the chain. In the scheme || such individual addressing 
is only required in the initialization. The method pro- 
posed in jj] is completely translation invariant in space. 
It requires homogeneous one- and two-local operations 
on 5-level systems. 

Here I describe a scheme for universal quantum com- 
putation via translation-invariant operations on a chain 
of qubits. No individual addressability is required. The 
scheme uses an Ising-type interaction and spatially uni- 
form one-qubit gates. The qubits are all of the same 
species. Cold atoms in optical lattices [|, where the re- 
quirement of local control adds to the overall technolog- 
ical challenge, are a candidate for the realization of the 
presented scheme. 



II. CONSTRUCTIVE ELEMENTS, 
UNIVERSALITY AND SCALABILITY 

Consider a one-dimensional chain of N qubits initial- 
ized in the state |00..0) which is repetitively updated ac- 
cording to the transition function 
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That is, in each elementary step of the evolution first 
a Hadamard gate is applied to each qubit and second, 
conditional phase gates are simultaneously applied to all 
pairs of neighboring qubits. This QCA transition func- 
tion has previously been discussed in |9j . 

Between the transitions one may apply translation- 
invariant unitary transformations of the form 



N 

U A {a) = (g)exp 

i=l 



(2) 



with A e {X,Y,Z}. (Note that the subscript l V labels 
the site. The same operation is applied to each qubit.) 
These requirements are equivalent to bang-bang pulses 
of form (J2J and a permanent Ising-type interaction 7i = 

Let us first observe that 
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where R is the reflection operator that sends the state 
of the qubit chain into its mirror image. Thus, despite 
the fact that the qubits at the end points need not be 
addressed, it is relevant that the chain has ends. (For 
adaption of the scheme to a ring of qutrits, see remark 1 
in Section ITTT|) . A proof of © is given in the appendix. 

Apart from its use in the computational scheme de- 
scribed here, the bit reversal operation R is interesting 
in its own. Recently, proposals for both approximate and 
perfect bit reversal in qubit chains with a Heisenberg- and 
XY- interaction have been made; See [I(J]-[l2| an d ref- 
erences therein. For perfect mirror reflection in an XY- 
chain the coupling strength needs to vary with position 
H3, but only mildly [3. 

Relation J3J represents a method to achieve spatial re- 
flection in systems with an Ising-type interaction. Here, 
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FIG. 1: The evolution of a local Pauli observable, Z3, for a 
qubit chain of length N — 8. The color-coded boxes denote 
Pauli operators Xi and Zi, respectively, and each row of boxes 
represents a tensor product of such Pauli operators. Within 
the cycle that leads to reflection of the chain, each local ob- 
servable undergoes the phases of expansion, transmission and 
contraction. When expanding or contracting, the operators 
pick up a sign factor of -1 under conjugation by a Y-pulse. 



the interaction strength is independent of position, but 
additional stroboscopic pulses (J2J) are required to realize 
the uniform Hadamard transformations. 

a. Universal set of gates. The key to the construc- 
tion of a universal set of gates for the described automa- 
ton is displayed in Fig. The N + 1-fold repetition 
of the elementary transition function constitutes half a 
clock cycle of the automaton. Within this period each 
local Pauli observable goes through the phases of expan- 
sion, transmission, and contraction to the mirror image 
of the initial position. During expansion and contrac- 
tion, the propagated observable is susceptible to a global 
V-pulse 



JY 



Y;=U Y (n)=Q<)Y i . 
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Namely, it picks up a sign factor under conjugation. Con- 
trarily, during the phase of transmission where the ob- 
servables behave as left- or right-movers [i| a F-pulse 
has no effect. The duration of expansion, transmission 
and contraction phase depend on the initial position of 
the local Pauli observable. Therefore, each local Pauli 
observable shows a characteristic response to sequences 
of Y"-pulses within the half-cycle. In this way, temporal 
control can be translated into spatial control. In the con- 
struction described below, suitably tailored sign-flips of 
Pauli observables are used to reverse rotation angles. De- 
pending on whether a rotation angle is reversed or not, 
two matching rotations will either cancel or amplify each 
others effect. 

The degree of spatial control obtained suffices to sim- 



ulate a quantum logic network. One constraint arises: 
every operation applied to the qubit chain is reflection- 
symmetric. Therefore, qubit i cannot be addressed sepa- 
rately from its mirror image at location i := N +1 — i. To 
cope with this constraint, the same network is simulated 
twice on the chain, once on the left side and once — as 
mirror image — on the right. This doubles the required 
length of the chain and also influences the readout pro- 
cess, as will be discussed in Section ITTlfr. 

To probe the available spatial control, consider the se- 
quence 

U z (c,a) = {TY CN TY CN - 1 T..TY Cl TY c °U z {-a/2)) 
(TY CN TY Cn - 1 T..TY C1 TY C ° U z (a/2) ) , 

(5) 

where c = (co,..,cjv) t is a binary vector. To analyze 
the effect of this sequence, the Pauli operators Y Ci are 
propagated backwards in time until they reach the Le- 
gates. There they accumulate to Y, 



Y = ^p— N Y CN T N ) (T~ 1 Y Cl T)Y' 
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where we have used that RYR = Y, and 



0, if [Y,Z l ] =0 

1, if {Y,Zi} = 



(6) 
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In this way, temporal control has been converted into 
spatial control, provided that — for suitable choices of c — 
the binary variables Si(c) do indeed vary with i. 

The Si are easily computed in the stabilizer formalism 
[13. Following [lJl, we write Pauli operators A in 

the form l) e r a , where a = ( ^ is a 2/V-component 

binary vector, z = («x, .., vn) t , x = (u>i, .., wn) T ', S, e £ 
Z 2 , and r a = <g)? =1 Z^Xp. The evolution of A under 
conjugation by our Clifford unitary T, A — > T(A) = 
TAT', may then be followed in terms of 6, e and a. 
The scalars 5 and e have no influence on the sign fac- 
tor (— l) Si JSJ, and we need to consider the update of a 
only, a(i) — > a(i+ 1) = Ca(t). Therein, C is a 2N x 27V 
binary symplectic matrix which takes the form 



C = 



r 1 
1 



(9) 



r is the adjacency matrix of the interaction graph (the 
line graph). Further, denote by F the 2N x 2/V-matrix 

I 

1 



F 



and observe that FC 



CF. Now, the 



vector s = (si, .., sn) t carrying the information about 
the sign flips under conjugation by Y is related to the 
vector c describing the temporal sequence of K-pulses, 
via s = Mzc. The matrix Mz encodes how temporal 
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control is converted to spatial control. Its elements are 
given by 



M z (i,t) 



a|*C*ay. 



(10) 



The binary quantities Mz(i,t) can be defined for all t 
(i = 1..N). But the matrix Mz is the collection of 
M z (i,t) only within the half-cycle < t < N. 

For the interval t £ [—1,N] (including [0,N]) the 
Mz(i,t) take the values 



M z (i,t) = 6(i — t — 1) + 9(t + i — N — 1). 



(11) 



Therein, the addition is modulo 2 and the step function 
is equal to one for all non-negative arguments, and zero 
otherwise. 

To prove one may — in addition to Mz — define a 
matrix Mx that encodes the sign factor acquired by the 
t-stcps-backward-propagated observable Xi under conju- 
gation with Y, Mx(i,t) :— a^.C'ay. Then, because 



of TX{T- X = Zi and TZ^ 1 = X % ® 
Mz(i,t) and M x (i,t) obey the recursion relations 



the 



M z (i,t+1) = M z (i,t-1)+ V M z (j,t), 

i|r„=i (12) 

M x (i,t+1) = M z (i,t). 

The solution of the recursion relation for Mz(i,t) is 
unique once the boundary conditions on two consecu- 
tive time slices are specified. In the discussed case, the 
boundary conditions are Mz(i, — 1) = Mz(i, 0) = 1, for 
all i = 1..N. The expression on the r.h.s of obeys 
the recursion relation (|12[1 in the interval < t < N — 1 
and the boundary conditions, and is thus the correct ex- 
pression for Mz- 

To perform a z-rotation on qubit i, one may apply a Im- 
pulse at times i — 1 and i. For ct = 5(t,i — l)+5(t,i), with 
((TTJl . one obtains Sj — S(j, i) + 5(j, N+l — i). Therefore, 
the pulse for a z-rotation of qubit i is 



= (T N+1 - i YTYT l - 1 U z (~a/2)) 
(T N+1 - i YTYT l - 1 U z {a/2)) , 



(13) 



with i = N + l — i the mirror site of i. In this way, 
two out of N z-rotations can be selected. A further 
discrimination between two sites i and i is not possi- 
ble because of the reflection symmetry of every sequence 
U(c, a) = RU(c, As a result of this symmetry, 

each quantum algorithm is run on the qubit chain in two 
copies, one being the mirror image of the other. This 
accounts for a factor of two in spatial overhead. 

From i|13|) the remaining universal gates can be de- 
duced easily. Relation (p~Sf) is conjugated by T _1 . Using 
T- X Z,T = X,, one finds 



e 2 2 



(T N ~ 

{T N- 



l YTYT l U x {-a/2)) 
i YTYT i Ux (ex /2)) . 



(14) 



a) 

physical 1 2 3 4 5 6 7 8 9 10 11 12 13 14 
position 

logical [1] [2] [3] [0][0] [3] [2] [1] 

•o#o#o#«o«o«o# 



I I mirror image 
qubits required in readout 



U z plZ) Y Y 



Z -rotation on qubit [2] 



V x (aJ2) Y Y 



X-rotation on qubit [2] 



U x (c</2) Y Y 



X¥-gate on qubits [1], [2] 



FIG. 2: Arrangement of logical qubits within the chain (a) 
and pulse sequences for gates (b), for n — 3. a) Only every 
second physical qubit in the chain, is used as a logical 
qubit. The qubits in between, "o", are ancillas. b) Pulse se- 
quences for the implementation of the universal gates, during 
the first half of a clock cycle. The second half-cycle is anal- 
ogous, with a — > —a. The symbol "■" denotes application 
of the elementary transition function T. The vertical bars 
denote pulses 



Conjugating relation (|14(l again by T 1 and noting that 
T- X X,T = Z, (g)^ r . , =1 Xj -. K, one obtains 



= T 



YTYT l U x {-a/2)T) 
T N - 1 - l YTYT l U x (a/2)T) . 



-iN-l-i 

-rjV-1- 



(15) 



Now, the length of the chain is doubled a second time. 
Specifically, the logical qubits in state \ip) are interlaced 
by ancilla qubits which remain in the state |0) throughout 
the computation, 

|^}l,2,3,..,n ► \lp') = 1^)1,3,5, ..,2n-l ® |0..0) 2 ,4,..,2n-2- 

(16) 

At this point it is suitable to introduce a logical coordi- 
nate [j] that is related to the physical location within the 
qubit chain via [j] = 2j — 1, for 1 < j < n, such that e.g. 

= V')[i],[2],..,[n] ® |0. .0)2,4, .., 2 „- 2 - Then, the action of 
the K%j gate is equivalent to a two-qubit next-neighbor 
entangling gate, 



exp 



(»'f W) = ex P Of x iJ] x [j+n) W)- ( 17 ) 



The gates (fH?)) , l(TI|) and (|T31) form a gate set that can be 
easily converted into the standard universal set [l6j |. 
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b. Readout. As an example for global readout I con- 
sider a measurement of the ^-component of the total spin, 
i.e., of the observable 

S z =Y,Zi. (18) 

i 

This is a model for readout e.g. of atomic qubits via res- 
onance fluorescence spectroscopy, if the atoms are well 
within the Lamb-Dicke limit (their separation and fluctu- 
ation of position are much smaller than the optical wave 
length) . In this case the underlying physics prohibits the 
extraction of local information. The readout method de- 
scribed below equally works for a scenario where a local 
readout could in principle be performed but is not pur- 
sued due to technological limitations. For specification of 
such a measurement model, see remark 3 in Section IlIII 

With the capability to perform arbitrary unitary evo- 
lution on a quantum register the total spin observable 
Sz acting on that register may be conjugated into any 
desired observable, and the standard network readout of 
the individual Zi should be feasible. However, for the 
setting described here the readout procedure is, in the 
case of probabilistic algorithms, slightly complicated by 
the fact that two copies of the algorithm are run simul- 
taneously and the readout measurement couples them. 

For the readout each of the two circuits on the chain 
requires one additional logical qubit that is in the state 
0) until the readout starts. The location of this qubit 
within the logical quantum register is denoted as [0] 
which shall correspond to the physical position 2n + 1 
in the qubit chain (see Fig. |2K)- The readout con- 
sists of three steps. First, Sz is measured and the 
outcome m is recorded. Second, for all logical qubits 
j = l..n the following procedure is performed: 2.1) 
Apply CNOT gates A(X) bm <8> A(X) wm ([0], [0] are 
the target qubits). 2.2) Measure the observable Sz 
and record the outcome m(j). 2.3) Apply the CNOT 
gates A(A)[ J ] [ ] (g> A(A)|^j^j again. Third, denote by 

J = {ji, .., those j < n for which m — m(j) = 2. 
If |J| < 2 the readout is finished. Otherwise, for all 
k = 2..\J\ perform the following procedure: 3.1) Apply 
the Toffoli-gates A 2 (X)y llUkm ®A 2 (X) I - ] J - m ([0], [0] 
are the target qubits), 3.2) Measure Sz and record the 
measurement outcome m(ji,jk)- 3.3) Apply the Toffoli- 
gates A 2 (X) [jl]t[jh] . m <g> A 2 (X) WhWm again. 

The conditional phase- and Toffoli gates are their own 
inverse such that the above protocol amounts to a se- 
quence of measurements where the first one is of Sz and 
the following ones are of conjugated observables. These 
observables mutually commute because they are all diag- 
onal in the computational basis, and the final state is a 
simultaneous eigenstate of them. The observables mea- 
sured in step 2 are Z[o]-%] + ^W%T+X^[o],|o] Z i- Since 
the qubits at locations [0] and [0] are individually in the 
state |0), and the observable Sz = J2i Z i already has 
been measured in step 1, step 2 effectively acts as the 
measurement of the observables + Zjjr. 



Similarly, step 3 amounts to the measurement of the 
observables H)^],^] (ll| + |ll)|-y|-y (11|, Vfc. This mea- 
surement assigns a sharp value to HiilHifc] "^HidHi*] - ^ s 
purpose is, in combination with the information gath- 
ered in step 2, to discriminate between the two cases 

i r Ui] = HifchHn] = Hi*]} and Hi] ^ Hi^l-Hii] ^r [jk] }, 
for each k. 

The state of the two quantum registers after the mea- 
surement, leaving out all the ancillas in the state |0), may 
be written as a superposition ^ Q c Q ,|r ( - a - > ) |f of compu- 
tational basis states |r( Q ))|r ( - Q - ) }, where r = (i"m, .., ?"[n]) T j 
r = (f]7p ■■' r Jn]) T - The measurements then impose con- 
straints upon which basis states are admissible, i.e., occur 
with nonzero coefficients c a . The constraints imposed by 
the measurements in step 2 read 

rj = f j = 0, for m ~ m(j) = 0, 

Hi] + Hi] m °d 2 = 1, for m — m(j) = 2, (19) 

Hi] = Hi] = 1 ' for in m(j) = 4, 

for all j = 1..N. The constraints that arise from the 
measurements in step 3 are 

Hiil+Hifc] mod 2 = 0, for m - m(j'i,ifc) = 2, 
Hii] + Hifc] mod 2 = 1) for m - m(Ji,jk) =0, 

for [ji], [j k ] £J,k = 2..\J\. The system ((TH1 EIB of linear 
equations has a unique solution if J = and two solutions 
otherwise. In the latter case, the two solutions r^ and 
r ( 2 ) correspond to two states \r^) \r^) and |r^ 2 ^ ) Ir^ 1 ^ ) , 
and the state of the quantum register after readout is 
their reflection-symmetric linear combination. The so- 
lutions r' 1 ) and r^ 2 ) are, modulo interchange, uniquely 
specified by the measurement outcomes m, {m(j)} and 
{m(ji, jk)}- For a probabilistic algorithm, the automa- 
ton produces two potential solutions in each run. 

c. Overhead. For the described method of univer- 
sal computation via translation-invariant interaction it 
is relevant that the elementary transition function T re- 
peated a sufficient number of times reduces to the iden- 
tity operation. Such an elementary 'clock cycle' takes 
2(N + 1) time steps where N is the length of the chain. 
But as described so far, one still ends up with a set of 
local and next-neighbor gates and thus pays the price of 
next-neighbor slowdown twice. This is not necessary. Be- 
fore the resources are counted, a long-distance entangling 
gate is described. As compared to a logically equivalent 
composition of next-neighbor and local gates, it cuts the 
required operational resources by a factor of order n. 

The natural gate set for the described scheme of quan- 
tum computation does not only contain local and next- 
neighbor gates, but all gates of the form 

U x (a, Lx,L 2 ) = exp (g) X^j . (21) 

As before, the logical position [I] corresponds to the phys- 
ical location 21 — 1 in the chain. The sequence for imple- 
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meriting Ux(oz, L±, L%) is 

U x {aMM) = (T t >YTYT t *U x (-a/2)T t i) . . 

(T ts YTYT t2 U x {a/2)T tl ) , [ ' 

withti = L 2 -Lx, t 2 = L x +L 2 -2 and t 3 = N -2Lx + 2. 
Now note that, with — ) := —X\—), 

e «r| >i!,j 2 < I = Ux(TT/2,h,l 2 )Ux(-TT/2,h + l,l 2 ) 
U x (7T/2,k + l,l 2 -l) 
U x {-n/2,h,l 2 -l). 

Thus, a long-distance x-controllcd spin-flip can be imple- 
mented with 4 elementary gates These four opera- 
tions can be grouped into two pairs, with constant value 
of L 2 in each pair. Either pair of operations can be imple- 
mented in one clock cycle because the F-pulse sequences 
of the two operations coincide. An x-controlled spin-flip 
thus takes two clock cycles irrespective of the distance 
between control and target qubit. 

Now counting the resources, a quantum computation 
on n qubits in the described scheme requires a qubit 
chain of length N = 4n + 2. Each clock cycle takes 
8n + 6 elementary time steps (applications of T), and 
each of the universal gates exp(«a/2Z/), exp (ia/2 Xi) 
and e i7r l ><i. ! 2\ I takes at most two clock cycles. As 
compared to the network model, there arises a constant 
overhead of 4 in the spatial and a linear overhead of 16n 
in temporal resources. 

III. CONCLUDING REMARKS 

Universal quantum computation can be performed by 
translation-invariant operations on a chain of qubits ini- 
tialized in the state |00..0). The described method re- 
quires translation-invariant interaction of Ising type and 
spatially uniform local rotations. As compared to the 
network model, there occurs a constant spatial and a lin- 
ear temporal overhead. 

Three remarks: 1) In the described setting trans- 
lation invariance is broken by the finite extension of 
the qubit chain. To create a setting with perfect 



translational invariance, one may replace the qubits 
by qutrits, where the additional level |2) docs not 
take part in any interaction, and consider a N + 1- 
qubit ring geometry instead of a chain. Set \^p{i)) — 
|0..0)i..j_i|2)j|0..0)i+i..jv+i. With either the translation- 
invariant superposition ^x+i ^h=i' I^W) or mixture 

W+l IvWX^K*)! as the initial state the scheme 

works as before. 

2) It may occur that a qubit chain capable of the re- 
quired interaction has been created but its length N is 
unknown. It can be found out easily by repetitions of 
the following protocol: Initialize the chain in the state 
|00..0), apply the transition function T t times and sub- 
sequently measure the spin observable Sz = $2»=i ^i- 
The received signal (Sz(t)} carries a characteristic signa- 
ture of AT, namely (S z (t)) = NS{t mod N + 1). 

3) Deviating from the setting discussed in Sectionlrflb. the 
available measurement may be such that local informa- 
tion could in principle be retrieved but is not retrieved, 
due to technological limitations. As an example, consider 
cold atoms in an optical lattice which are read out via 
resonance fluorescence spectroscopy and whose separa- 
tion is larger than the wave length of the probing laser. 
In this case, the readout can be modeled by the map 

p^p'(m)~ Yl P({m i })pP({m i }), (24) 

{ m dl m i~ m 

with the projectors P({nii}) = ® i=1 1+ ™ iZi , and mi = 
±lVi = 1..N. Every atom is individually measured but 
only the global signal m = J^- m, is recorded. For mea- 
surements of this type the readout procedure is the same 
as discussed in Section |H] and has the same efficiency. 
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and a relation that expresses x in terms of z. From 
taking into account that T represents a line graph, 
the following recursion relations are obtained for z = 
(vi, ..,v N ) T and x = (wi, ..,w N ) T : 



«i(t+l) = «a(t)+t;i(t-l), 

Wi(* + 1) = (t)+«<+i (*)+«<(* 

V N (t + l) = V N -t(t) + VN{t - 1), 



1), 1 < i < N, 



APPENDIX A: BIT REVERSAL 

Sufficiently many applications of the elementary tran- 
sition function T reverse the order of qubits within the 
chain. As stated in ®, with R the spatial reflection, 

T N+1 = R 

This phenomenon is related to the transmission of one- 
qubit states through chains described in ^tJ- Rela- 
tion 121 is now proved. What needs to be shown is 
T N+1 (X p ) = X ¥ and T N+1 (Z p ) = %, for all p = 
1, . . , N . For this purpose, the vector space formulation 
of the stabilizer formalism is used. In particular, 
I use the conventions and results of 15]. The evolution 
of the Pauli observables X p , Z p is followed in terms of 
the binary quantities a and e, introduced in Section [n] 
below Eq. ©. 
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FIG. 3: The solution for the recursion relations (IA3II in free 
space. Each vertex in the lattice denotes a space-time point 
(i,i). The symbol "•" indicates Vi(t) — 1 and "+" Vi(t) = 0. 
To obtain the solution for the case with boundary from the 
free space solution, the lattice is 'folded back' to the space 
interval [0, N + 1]. 



First, the recursion relation a(t + 1) = Ca(f) is trans- 
lated into a recursion relation for the z-part of a alone 



Wi(t+1) = Vi{t), l<i<N. 

(Al) 

We seek the solution of these recursion relations for the 
time t = N + 1 , with boundary conditions 



Vi{-1) = 0, Vi(0) = 5(p-i), VI < i < N. 



(A2) 



The translation-invariance of the recursion relations (|Alj) 
is broken by the finite extension of the chain. The prob- 
lem of solving these recursion relations is now reduced to 
the translation-invariant case of the infinite chain. Note 
that if a configuration {£>;(£)} obeys the recursion rela- 
tions 



Vi(t + 1) = Vi-i(t) + Wj+i(t) +«j(t - 1), 



then 



Vi(t) 



s<+2i(iv+i) 

/ — — oo 



-i+2/(jV+l) 



(A3) 



(A4) 



obeys the recursion relation (|A1|) . The reduction of the 
case with spatial boundary to the case without boundary 
is illustrated in Fig. |3 
It is easily checked that 

Vi(t) = {9(p - i + 1) + 9(r - i - 1 - 1)) 5(p-i+t mod 2) 

(A5) 

is the solution of 1A3(1 with the boundary conditions 
Vi{— 1) = 0, Vi(0) = 5(p — i), Vi. If this expression is 
inserted into the r.h.s. of (| A4|> one obtains a solution for 
Vi(t) with boundary conditions (|A2|) . It is observed that 



Vi{t) 
Wi(t) 



6( P 
S(p- 



i + t mod 2), 
i + t + l mod 2). 



(A6) 



Next, the solution v(t) 
N + l. One obtains 



5j is evaluated for t 



Vi(N+l) = S(p-(N + l-i)). 



(A7) 



This implies T N+1 (Z P ) 

Next, the sign factor (— l) e in the above relation is 
worked out. Using the result of [l5|, the recursion rela- 
tion for e reads e(t + l) = e(t) + z(t) T T L z(t) + x{t) T z(t), 
with Tl the lower triangular part of T. With l|A6jl . 



C(i) = 6(0) 

with T(X p ) 



such that T N+1 (Z p ) = Z p . Finally, 



Z p , T N +\X p ) 



T^T N+1 T(X p ) 



X p . 



